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ABSTRACT 

We suggest that the deceleration of the relativistic shock by a denser part of 
the interstellar medium off line-of-sight produced the observed radio "flare" in 
the early afterglow of GRB990123. We find that this scenario is consistent with 
observations if the particle number density of this denser part of the medium is 
between ~ 200 and ~ 2 x 10 4 cm" 3 . Because of the premature deceleration of 
part of the shock, the later stage of the afterglow should decay modestly faster 
than the powerlaw expected from an isotropic shock propagation. 

Subject headings: gamma rays: bursts — cosmology: miscellaneous 
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1. Introduction 

GRB990123 had some remarkable characteristics that drew intense interest from 
astrophysicists. Its source appears to be at a cosmological redshift z > 1.61 (Kelson et 
al. 1999). Assuming isotropic emission this leads to a huge energy release: ^ 3 x 10 54 
ergs in 7-rays alone (Kulkarni et al. 1999a). Prompt optical follow-ups saw a very bright 
(magnitude 9) early afterglow optical flash (Akerlof et al. 1999). Since the 7-ray burst 
(GRB) itself, its afterglow has been monitored almost constantly in X-rays, the optical 
band and radio. 

Radio observations of GRB990123 have been particularly puzzling. As expected from 
afterglow models, initial measurements (six hours after the burst) obtained only upper 
limits (Frail et al. 1999a). One day later, however, new observations showed a 8-cr detection 
at 8.46 GHz of 260±32/iJy (Frail et al. 1999b). On the other hand, observations at 4.88 
GHz made during a period overlapping with the 8.46 GHz detection, yielded only upper 
limits (Galama et al. 1999). Finally, three to five days after the burst the afterglow radio 
output was again consistent with zero (Kulkarni et al. 1999b). Figure 1 summarizes the 
detection and upper limits. The radio emission is rather peculiar, as it does not conform to 
the gradual rise-and-decay time profile with a timescale of O(10) days which is expected 
from the standard GRB afterglow model and which was rather successfully confirmed by 
the radio afterglow observations of GRB970508 (Waxman, Kulkarni and Frail 1998). 

Several possibility have been raised to explain this one-time radio "flare" . Kulkarni et 
al. (1999b) have suggested that the flare is an interstellar scattering and scintillation event. 
Alternatively, it might be part of the early afterglow from the external reverse shock, if 
self-absorption of the radio emission suppressed the radio flux during the first day as well 
as during the second day but only at the 4.88 GHz band (Sari and Piran 1999). 

Here we investigate a third possibility (Shi and Gyuk 1999), that the radio "flare" 
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might be due to the relativistic shock ploughing into a dense part (a cloud, or ejecta, for 
example) of the interstellar medium (ISM) off line-of-sight (LOS). This off LOS portion 
of the shock was therefore decelerated more efficiently than the rest (including that along 
LOS), and in turn gave rise to the premature emission of radio which also faded away 
relatively rapidly. We envision a geometry as in Figure 2. We will refer generally to this 
denser part of ISM off the LOS as a "cloud" . 

In the standard GRB afterglow model, the afterglow is generated by synchrotron 
radiation in the external shock of the GRB event as the shock gradually decelerates in a 
homogeneous ISM (with density n ~ 1 cm -3 ). The frequency of the synchrotron radiation 
depends strongly on the Lorentz factor, j e , of the electrons in the shock, which in turn 
scales with the Lorentz factor, 7, of the shock. Therefore, the afterglow starts at shorter 
wavelengths, in X-rays, and progressively moves to longer wavelengths, as the shock 
gradually slows down in the ISM. After O(10) days, the afterglow peaks in the radio band. 
A similar progression is envisioned to occur in the portion of the shock that encounters 
the off-LOS cloud, except over a much shorter timescale. However, because the afterglow 
radiation is beamed to an opening angle ~ I/7 we can only see the off-LOS cloud shock if it 
is within ~ I/7 of the LOS. Thus relativistic beaming alone may prevent short wavelength 
(optical etc.) "flares" originating from the off-LOS cloud from being detectable. 

If the size of the cloud is comparable to the distance of the cloud to the GRB source, 
we can crudely approximate the deceleration of the relativistic shock in the cloud as if the 
shock were decelerating in an homogeneous ISM of enhanced particle density {n ^> 1 cm~ 3 ) 
and with spherical symmetry. This approximation should hold sufficiently well for the later 
epoch of the deceleration, which is relevant to radio emission. In so doing, we employ the 
scaling relations in the standard afterglow model. 
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The Lorentz factor of the relativistic shock scales as 



J PS 




(1) 



from energy and momentum conservation considerations (see e.g., Piran 1998). In equation 
(1), E 52 is the initial energy of the shock in units of 10 52 erg, assuming a 4tt expansion 
angle; n\ is the particle number density of the medium in cm -3 ; 7100 is the initial Lorentz 
factor of the shock in units of 100; tday is the time elapsed since the GRB in days, as seen 
by the observer; and z is the redshift of the GRB. In the radiative regime, the particles 
in the shock convert their kinetic energy into radiation rather efficiently In the adiabatic 
regime, the radiation loss is negligible. 

For the external shock in a GRB, the transition from the radiative regime to the 
adiabatic regime occurs at a time (Piran 1998) 



where e e is the fraction of thermal energy of the shock that resides in the random motion of 
electrons, and e# is the ratio of the magnetic field energy to the thermal energy density of 
the shock (~ A^ 2 riim p c 2 where m p is the proton mass and c the speed of light). Canonical 
values are e e ~ 0.6 and e B ~ 0.01, obtained by fitting the standard afterglow model to 
the observed afterglow of GRB970508 (Wijers and Galama 1998; Granot, Piran and Sari 
1998). Before t r ^ a the cooling time is shorter than the dynamic timescale, and the shock is 
radiative. 

For an energetic GRB event such as GRB990123 where E52 ~ 10 3 , and a cloud much 
denser than the average ISM (rii 1), the transition occurs much later than a day. We 
therefore assume a radiative shock in our treatment. 

There are three synchrotron emission frequencies that are crucial: v m , the peak 
synchrotron radiation frequency if electrons in the shock are slow-cooling; v c , the peak 
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synchrotron radiation frequency if electrons are fast-cooling; and u a , the self- absorption 
frequency of the synchrotron radiation below which the radiation is absorbed by electrons 
in the shocks. Depending on the ratio of these key frequencies, the synchrotron radiation 
from the external shock can have very different spectral shapes. 

In the standard afterglow picture, the shock-heated electrons develop a power law 
number density distribution N(j e ) oc % p where 7 e > r y e ,min is the Lorentz factor of 
electrons. The minimum Lorentz factor cut-off is 

7e >mi n « ^^e7 ~ 2.2 x 1 3 (ejO^E^n^ Vo^dlfK 1 + *)/ 2 -6] 3/7 , ( 3 ) 

where m p and m e are proton and electron masses respectively (Piran 1998). The power law 
index p is found to be ~ 2.5 by fitting the observed GRB spectra and that of the afterglows. 
If electrons are slow-cooling, their peak synchrotron emission will be in the observer's frame 
at a frequency 

u m = 7 ^7 2 ^ eC « 7x 10 u (6 fl /0.01)V' (Ce ^ 

(4) 

where e is the electron charge and B is the magnetic field. 

If, however, the shocked electrons cool quickly, they will mostly radiate from a cooled 
state, whose Lorentz factor is 

7 - ~ 4^8^ ^ 7 X 104 (eB/OMrE^W^l tZ [(1 + z)/2.6]-^ (5) 

where ot = 87re 4 /3mgC 4 = 6.65 x 1CT 25 cm 2 is the Thompson scattering cross section (Piran 
1998). The emitting frequency in the observer's frame is 

+ 



T~~~7^~~ - 5 x 10 15 (es/Omr^E^n^^tSt^ 7 [(1 + ,)/2.6]- 5 / 7 Hz. (6) 
1 + z InmpC y 



To find the self-absorption frequency u a , we follow Granot et ed. (1999) and Wijers and 
Galama (1999) to calculate at what frequency the optical depth becomes unity. A crude 
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estimate of the optical depth r is r ~ a' u ,R/^, where a' v , is the absorption coefficient at a 
a frequency u', and R/j is the thickness of the shock, all in the local rest frame. However, 
we cannot directly adopt the formula for v a in Granot et al. (1999), or Wijers and Galama 
(1999), because both have assumed a slow electron cooling regime in an adiabatic shock. 
In our problem, the electrons are fast-cooling, and the shock is radiative. We therefore 
substitute the electron Lorentz factor 7 e c (fast-cooling regime) for 7e im j„ (slow-cooling 
regime), and likewise substitute the shock Lorentz factor 7 in radiative shocks for that in 
adiabatic shocks. The resultant self-absorption frequency is then 



We have implicitly assumed v a <^ v c since for rapid cooling most of the electrons will be 
at Lorentz factor 7 e c . For these electrons, absorption at frequencies 3> v c falls off rapidly. 
If we further assume that the time-integrated absorption from newly injected electrons in 
transition to their final cooled state is small, we will always have v a <; v c . 

Therefore, with E 52 ~ 10 3 , n\ >> 1 and tday ~ 1, and canonical values of e#, e e and 
7100, there is a hierarchy in frequencies: v a <; u c < u m . The peak flux of the synchrotron 
radiation seen by an observer is at u c (Piran 1998): 



FM) = F v , max w 1.8 x 10 3 .F^/O.Ol) 1 / 2 E% 7 4 lw 8 J 7 1^ 7 [(l + z)/2.6] 10 / 7 d 28 2 ^Jy, (8) 



where T < 1 is a geometric factor to account for the fact that the emission is from off LOS 
so we only see an edge of the radiation cone, and d 2 s is the luminosity distance in the units 
of 10 28 cm. Fluxes at other frequencies are (Piran 1998) 



10 8 (WO.OI) 6 / 5 E'£ ni 7K0 C 5 K 1 + ^/S.e]" 1 / 5 Hz. 



(7) 



F 




if v a < v < v c \ 



(9) 



nu 



{y/v a )-V*F v {v c ) 
{vlv a )-vl 2 F v {v m ) 



if v > v, 
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The spectrum and its evolution is schematically plotted in Figure 3. 

We require u c , v a > 8.46 GHz at the time of td ay ~ 1 so that a non-detection of 
radio signal at 4.88 GHz is compatible with a simultaneous detection at 8.46 GHz. This 
condition is satisfied if ri\ ^ 2 x 10 4 . Assuming at a redshift z ~ 1.61, the luminosity 
distance to GRB990123 is g?28 ~ 4 (its order of magnitude is insensitive to different 
choices of cosmology). The peak flux is then F^ max ~ 3 x 10 5 Fn^ 14 /xJy when adopting 
for other parameters values mentioned above. Scaling down to 8.46 GHz, we find 
F„(8.46 GHz) ~ F„ )mas (8.46 GHz/z/ c ) 2 ~ 3 x 10~ 3 Fnf' 14 /iJy.Q To yield a 260 //Jy detection 
would therefore require n\ ^ 200. Because at this part of the spectrum F v oc z/ 2 , a flux of 
260±32 /iJy at 8.46 GHz implies a flux of 86 /xJy at 4.88 GHz, consistent with the 3a limit 
of 130 //Jy measured at this frequency (Galama et al. 1999). Given 200 <^ rt\ ^ 2 x 10 4 , the 
Lorentz factor of the prematurely decelerated portion of the shock at t^ay ~ 1 is of order 
G(l). 

The non-detection of radio emission six hours after the burst may be due to strong 
absorption (i.e., u a too large), or it may simply be that the shock hadn't yet encountered 
the cloud. While three days later, this portion of the shock has become very weak, and its 
emission is further absorbed by the main shock that propagates along LOS. It should be 
kept in mind that a factor of several below the level of the detected emission might render 
the emission undetectable. 

Assuming that the dimension of the cloud is comparable to the size of the fireball 
~ 47 2 t at tday ~ 1 ; w e find a mass for the cloud to be ~ 10 -5 to 1O~ 3 M . We speculate 
that it may be ejecta from the GRB site. 

The main portion of the relativistic shock along LOS is not affected by the cloud 
l We have assumed v a ~ u c , which will be the case for n\ ^ 200. 
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off LOS. It generates the main afterglow as expected from the standard GRB afterglow 
model. The temporal decay of this afterglow in a given frequency band follows a powerlaw 
t L1 (Vietri 1997; Waxman 1997; Sari, Piran and Narayan 1998). But as its radiation 
cone become wider (opening angle 9 ~ I/7), the viewing area of an observer is larger. 
Eventually the area will engulf the portion of the shock that was prematurely decelerated 
and terminated. The temporal decay of the main afterglow should then be faster than 
t~ 1A . The transition to a faster decay law is not unique: for example, if the relativistic 
shock is a narrow jet, the temporal decay of its afterglow steepens when its opening angle 
9 < I/7. Depending on the details of the model it may either steepen by an additional 
£-3/4 p 0wer (Meszaros & Rees 1999), or steepen to t~ p (Kulkarni et al. 1999). The rate of 
afterglow decay due to an off-LOS hole in a spherical shock should be more modest than 
that due to a jet and indeed in this scenario, we should expect the afterglow decay will 
eventually approaches its initial shallower decay profile, as the influence of the geometric 
defect becomes increasingly less significant. 

In summary, we show that the radio "flare" observed in the early afterglow of 
GRB990123 may be due to a relativistic shock encountering a denser part of the ISM, (with 
a density between ~ 200 and ~2x 10 4 cm~ 3 ) off line-of-sight. A transition from a t~ 1A 
decay to a modestly faster temporal decay is expected in the later stage of the afterglow. 
This scenario also implies that the relativistic shock that generates the afterglow is unlikely 
to be beamed by more than a factor of a few. 

We thank George Fuller, Bob Gould and Art Wolfe for discussions. XS acknowledges 
support from NSF grant PHY98-00980 at UCSD. GG wishes to thank the Department of 
Energy for partial support under grant DEFG0390ER40546 and Research Corporation. 
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Fig. 1. — Radio observations at 8.46 GHz (the solid lines, upper limits are 2a, Kulkarni et 
al. 1999b) and 4.88 GHz (the dashed lines, 3a upper limits, Galama et al. 1999). 
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Fig. 2. — A possible geometry near the site of GRB990123. 
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Fig. 3. — A schematic plot of the afterglow spectrum from the off line-of-site deceleration in 
the cloud. 



